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One can confine the two-dimensional electron 
gas in semiconductor heterostructures electro- 
statically or by etching techniques such that a 
small electron island is formed. These man-made 
"artificial atoms" provide the experimental real- 
ization of a text-book example of many-particle 
physics: a finite number of quantum particles in 
a trap. Much effort was spent on making such 
"quantum dots" smaller and going from the meso- 
scopic to the quantum regime Far-reaching 
analogies to the physics of atoms, nuclei or metal 
clusters were obvious from the very beginning: 
The concepts of shell structure and Hund's rules 
were found to apply just as in real atoms! 
In this Letter, we report the discovery that elec- 
trons confined in ring-shaped quantum dots form 
rather rigid molecules with antiferromagnetic or- 
der in the ground state. This can be seen best 
from an analysis of the rotational and vibrational 
excitations. 

While the independent-particle picture was successful 
in describing the electronic structure for rather large par- 
ticle densities, for more dilute systems or in stronger mag- 
netic fields correlation effects are of crucial importance. 
Configuration-interaction (CI) calculations, which have a 
long tradition in quantum chemistry and cluster physics, 
were then much used |3|. Although these so-called "ex- 
act" calculations arc numerically demanding and limited 
to the smallest sizes, they still are able to provide signif- 
icant insight into the many-body phenomena that occur 
in these finite fermion systems with reduced dimensional- 
ity. In this Letter we apply CI techniques to investigate 
the electronic structure of quantum rings that contain 
up to seven electrons. Usually the confinement of small, 
two-dimensional quantum dots is to a very good approx- 
imation harmonic. Correspondingly, we model quantum 
rings as they are realized in the laboratory by a poten- 
tial of the form V{r) = ^m*ujQ{r — rp)^. For moderate 
confinement this potential corresponds to an harmonic 
dot with its center removed. Ground and excited states 
of N electrons trapped in the potential V{r) are deter- 
mined from numerical diagonalization as a function of 
the total angular momentum. Surprisingly, at electron 
densities and strengths of the ring confinement where one 
should expect electron liquid behavior, a model which as- 
sumes localization of the electrons in the ring is successful 
in analyzing the many-body spectra. Group-theoretical 
methods familiar from molecular physics provide the nec- 
essary tools to uncover rotational and vibrational struc- 
tures in the spectra. The spin sequence and energies of 
the low-lying states for given angular momentum can be 



understood from the symmetry associated with the elec- 
tronic ground state configuration. It is intriguing that 
the success of the simple rigid-rotor model for the low- 
lying states is not limited to a regime where the system 
becomes very one-dimensional. The fact that the elec- 
trons behave as if they were localized in the ring is also 
reflected in a remarkable agreement of the CI results with 
the Hubbard or Heisenberg model. Localization at large 
electron densities has earlier been discussed in parabolic 
quantum dots, where the interpretation is not yet con- 
clusive B. We write for the Hamiltonian 
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where m* and e are the effective mass and the dielec- 
tric constant of the corresponding semiconductor mate- 
rial. The parameters that determine the properties of the 
quantum ring are the number of electrons N , the radius 
of the ring tq and the strength wq of the harmonic con- 
finement in the radial direction. The quantities tq and 
Wo are related to the more fundamental quantities Ts, 
the one-dimensional density parameter which describes 
the particle density n = I/(2rs) along the ring (thus 
To = Nts/tt) and Cp, a dimensionless parameter that 
measures the degree of one-dimensionality. Cp essen- 
tially describes the excitation energy of the next radial 
mode hujo, which is defined to be Cp times the (ID) 
Fermi energy. We thus obtain hujo = Cph'^n'^ /{^2m*r'^J . 
The higher the value of Cp, the more the radial modes 
are frozen in their ground states. Thus, the ring is nar- 
rower for larger Cp. For the CI calculation, the spatial 
single-particle states of the Fock space are chosen to be 
eigenstates of the single-particle part of the Hamiltonian 
H. We expand them in the harmonic oscillator basis. Ac- 
cording to their eigenenergies, from 30 to about 50 lowest 
single-particle states are selected to span the Fock space. 
Typically this means that for lower angular momentum 
states several radial quantum numbers n = 0, 1, 2, 3 are 
included, whereas the higher angular momentum states 
/ = ±6, . . . , ±10 have only n — 0. To set up the Fock 
states for diagonalization, we sample over the full space 
with a fixed number of spin down and spin up electrons, 
Ni + — N . From this sampling, only those states 
with a given total orbital angular momentum and a con- 
figuration energy (corresponding to the sum of occupied 
single-particle energies) less than the specified cutoff en- 
ergy Cc are selected. The purpose was to choose only the 
most important Fock states from the full basis and hereby 
reducing the matrix dimension to a size d ~ 2-10^. To ob- 
tain all the eigenstates with different total spin, we have 
to set iVj = iV| = N/2 for even particle numbers (5*2 — 0, 
all states with different total spin have this component), 
and analogously iV| = ± 1 for odd numbers. Once 
the active Fock states have been specified, the Hamil- 
tonian matrix is calculated. For diagonalization we use 
the Arpack library ||^ suitable for large, sparse matrices. 
Finally, the total spin of each eigenvector is determined 
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by calculating the expectation value of the S'^ operator. 
The many-particle states are characterized by the total 
orbital angular momentum M and the total spin S. The 
lowest energy at given angular momentum M defines the 
so-called "yrast" line, a terminology that was introduced 
in nuclear physics many years ago and comes from the 
Swedish word for "the most dizzy". 

To begin with, let us look at the energy spectra of 
a ring with six electrons at a density corresponding to 
= 2ag. (Throughout this paper, we use effective 
atomic units. Taking GaAs as an example, the units of 
length and energy are ag = 9.8nm and a.u.* = 12meV.) 
For Cf = 4, the strength of the ring confinement is mod- 
erate, as it can be seen from the density distribution dis- 
played in Fig. 1^. 
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FIG. 1. Density along the radial direction for A'^ = 6 and dif- 
ferent parameters rs and Cf- The confinement is moderate for 
Vs — 2 and Cf = 4 {green). The degree of one-dimensionality is 
enhanced for Cf = 10 (red). The low-density regime is displayed 

for Ts = 6 and Cf = 25 (blue) 
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FIG. 2. Many-body spectra of a quantum ring with A'^ = 6 
electrons (r^ — 2ag and Cf = 4). The energy difference be- 
tween the yrast states for AI — and M = 6 is 0.143 a.u.*. 
The spin S is given for the low-lying states (as also indicated by 
different colors of the levels) 
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The system parameters TV, Vs and Cf in this case are 
chosen such that the rings resemble quite closely struc- 
tures as they can nowadays be made in the laboratory. 
For the six electron ring, Fig. ^ shows the 50 lowest states 
for all angular momenta from M = up to Af = 6. The 
spin configurations are given for the low-lying states. The 
ground state for M = has spin S — and is followed 
by a state with 5 = 2, then S — 1 and, close in energy, 
another S ~ state. A large energy gap separates these 
low-lying states from higher bands with a much increased 
density of states. We note that the sequence of spins of 
the states below the gap for M = is repeated at M = 6. 
(The energy difference between the two 5 = states at 
M = 6 is slightly reduced, as in a not very narrow ring ro- 
tation expands the ring.) Inspecting the grouping of the 
states below the gap for each of the different M-values 
more closely, we see that in a similar way, the states for 
M = 1 follow the same sequence as those for M = 5, and 
the M — 2 spectrum is repeated at M = 4. 



FIG. 3. Spectra of a narrow low-density quantum ring with 
six electrons [rs = 6 and Cf = 25), showing the seven lowest 
vibrational bands. For the vibrational ground state we show the 
M = levels in an magnified scale to demonstrate that the en- 
ergy ratios are the same as in the much wider rings shown in 
Figs. 2 and 5. 

These sequences of states clearly reveal the spectrum 
of a dynamical system that can be described, at least 
approximately, in terms of independent rotation, vibra- 
tion and intrinsic spin degrees of freedom. This becomes 
particularly clear when repeating the above calculation 
for a much lower particle density in a more narrow ring. 
We show in Fig. || the 6-particle spectra for = 6a*g 
and Cf=25. (The particle density for these parameters 
is displayed in Fig. |^). The spectra now indeed consist of 
narrow bands. The lowest band is the vibrational ground 
state and consist of the rotational levels. The different 
spin levels become almost degenerate, as localization has 
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reached a degree where spin-spin interactions become less 
important. Nevertheless, the sequence of spins and their 
relative energetic order is similar to the higher density 
results shown in Fig. (1), as a plot of the low- lying levels 
at an enlarged energy scale for M = demonstrates. In 
the following we analyze the above spectra in terms of 
the standard effective hamiltonian employed in the in- 
terpretation of the rotational and vibrational spectra of 
planar polygonal molecules composed of N identical spin 
1/2 fermionic atoms |5|; the standard description is very 
slightly modified in order enforce the planarity and the 
one-dimensionality of the molecule considered here: 

N 

Hef f = AM^ + ^^aUa + J S, • S, (2) 

where A is the rigid moment of inertia of the N " atoms" 
located at the N vertices of an equilateral polygon, uja 
are the frequencies and Ha (= 0, 1,2,.. .) the number of 
excitation quanta of the different normal modes of vibra- 
tions; for odd N there are ^(A^ — 1) two fold degenerate 
normal modes, while for even N there are ^A^ — 1 twofold 
degenerate and one non-degenerate modes. The last term 
in Eq. (^) describes the nearest neighbor spin-spin inter- 
action that is a result of the exchange term in the inter- 
action of the neighboring fermions. In the expression (|^) 
we have ignored the contribution of a possible spin orbit 
term. 

For the N = 6 molecule discussed in Figs. || and ||, the 
Cqv symmetry classification of the normal modes is (in or- 
der of increasing vibrational frequency) Ei (2-fold degen- 
erate), i?2 (2 fold degenerate) and Bi (non-degenerate). 
The rotational states carry symmetry for M = {Ai or 
A2), M EE 1 or 5 {El) M = 2 or 4 {E2), M = 3 (Bi or 
B2), where the congruent sign in all these terms refers 
to congruent (mod 6), and the total wave function must 
have Cqv symmetry B2 in order to fulfill the requirement 
of anti-symmetry between all the electrons. 

As we learned from the above analysis, the lowest band 
is of rotational nature, while the bands above the en- 
ergy gap represent vibration in addition to rotation. The 
three vibrational modes of such a system have energy ra- 
tios wi : 0^2 : 1^3 = 1 : \/7/3 : \/3, in excellent agreement 
with the three lowest excitations for M = shown in 
Fig. H The three highest states have then energies 2uji, 
LUi + UJ2 and LUi + lus- 

A generic model for localized electrons is the Hubbard 
model 1^,^ which can be solved exactly for a small num- 
ber of lattice points. The real electrons in the ring inter- 
act with the long-range repulsive Coulomb force. Clas- 
sically the electrons would form a ring of equidistantly 
localized electrons. In the Hubbard model this case corre- 
sponds to the so-called half-filled model, where the num- 
ber of electrons equals the number of lattice points. In 
the strong interaction limit it can be transformed to the 
anti-ferromagnetic Heisenberg model |@] leading to the 
last term of Eq. (^. 
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FIG. 4. Schematic view of the spin configurations in a 
six-electron ring for different values of S. 



On the basis of the results from exact diagonalizations 
and their group-theoretical analysis described above it is 
now interesting to see how well such simple models of 
localized electrons can describe the many-body spectra 
in these finite, not truly one-dimensional quantum rings. 
Fig. ^ compares the energy spectra of the Heisenberg 
model to the CI spectra described above. In the latter 
case the center-of-mass motion associated to the orbital 
angular momentum (AM^-term) is subtracted from each 
level. The coefficient A is determined to make the low- 
est energy eigenvalue corresponding to Af = equal to 
that of M = 0. The resulting values of A are very close 
to A — TT^ / (2N^r1) obtained for strictly localized elec- 
trons. The scale of the Heisenberg result is determined to 
get the same band width as in the case of real electrons. 
Fig. 1^ displays for electron numbers = 5, 6 and 7 at 
rg — 2a|j and Cf = 10 the spectra of the states below the 
first gap which (as analyzed above) all are of rotational 
nature. For N = 7, corresponding to the larger particle 
number, there are now more states as we have more pos- 
sibilities for the different spin configurations. The above 
comparison shows that indeed, the electrons in a ring can 
be rather accurately described as localized electrons: The 
whole rotational spectrum close to the yrast line can be 
determined by a spin model combined with a rigid center- 
of-mass rotation. 

It should be noted that the results shown in Fig. ^ 
are calculated for a ring with a finite width. Making 
the ring narrower the small disagreement between the 
Heisenberg model and exact calculation vanishes. We fi- 
nally mention that the pair correlation function does in- 
dicate anti-ferromagnetic coupling of the electrons in the 
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ground state. (See configuration schematically shown in 
the top left corner of Fig. ||). However, the correlation 
is not as clear as expected from the rotational structure 
in the energy spectrum. In the limit of an infinitely long 
one-dimensional Heisenbcrg model the spin-spin correla- 
tion decreases as l/j'ij. Consequently, even in the more 
narrow rings, where the rotational and vibrational spec- 
tra show a clear localization of electrons, the pair corre- 
lation function shows only rather weak spin-spin correla- 
tion. 




N=5 



— 3 — 



3 — 


7/2 


2 — 


5/2 


1 — 


3/2 


— 


1/2 



— Hubbard 



M=0 —Angular Momentum— ^ M=5 



FIG. 5. Low-lying (rotational) spectra of quantum rings with 
A'^ = 5, 6 and 7 electrons: Comparsion between the Heisenberg 
model (left columns in each of the spectra for given M) and 
exact diagonalization (right). The spins of the levels are marked 
by different colors 

In conclusion, we reported rotational and vibrational 
many-body spectra of quantum rings confining up to 
seven electrons. We saw that the results obtained by 
configuration interaction calculations are best described 
by assuming localization in the internal structure of the 
ground-state many-body wave function. For even elec- 
tron numbers, antiferromagnetic ordering was found for 
the ground state. Such so-called spin density waves in 



quantum rings actually were predicted by density func- 
tional theory The many-body spectra of the re- 
ported exact diagonalization calculations, their analysis 
by group-theoretical methods and further the comparison 
to the Heisenberg model confirm that indeed the simple 
mean-field picture (as it is provided by density functional 
theory) can to a rather large degree correctly map out the 
internal symmetry of the many-body wave function. 
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